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■ ABSTRACT 

^< \ Some formulae are presented for finding two-integral distribution functions (DFs) 

■ which depends only on the two classical integrals of the energy and the magnitude 
OA | of the angular momentum with respect to the axis of symmetry for stellar systems 

with known axisymmetric densities. They come from an combination of the ideas of 
\ Eddington and Fricke and they are also an extension of those shown by Jiang and 

Oh Ossipkov for finding anisotropic DFs for spherical galaxies. The density of the system 

^ (J-i is required to be expressed as a sum of products of functions of the potential and of the 

radial coordinate. The solution corresponding to this type of density is in turn a sum 
of products of functions of the energy and of the magnitude of the angular momentum 
about the axis of symmetry. The product of the density and its radial velocity disper- 
sion can be also expressed as a sum of products of functions of the potential and of the 
radial coordinate. It can be further known that the density multipied by its rotational 
velocity dispersion is equal to a sum of products of functions of the potential and of 
the radial coordinate minus the product of the density and the square of its mean 
rotational velocity. These formulae can be applied to the Binney and the Lyndcn-Bcll 
• models. An infinity of the odd DFs for the Binney model can be also found under the 

| assumption of the laws of the rotational velocity. 
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1 INTRODUCTION 



?—j ' It is the most straightforward to construct self-consistent stellar systems by means of finding distribution functions for a 
. . i stellar system with a known gravitational potential. Once the potential of the system is known, the mass density p of the 
system can be uniquely detemined via Poisson's equation generated by the well-known Newtonian gravitational law, and the 
structure of the stellar orbits can be also obtained according to Newton's equations of motion. The system is hence constructed 
from the structure of the orbits in the potential. This construction is also the so-called "from p to /" approach for finding a 
self-consistent distribution function / (Binney & Tremaine 1987, hereafter BT; Hunter & Qian 1993). Since the mass density 
is the integration of the distribution function (hereafter DF) over the velocity variable in the phase space of the system, the 
problem of finding the DF is that of solving an integral equation. In a system with a known spherical potential, there is a 
unique isotropic DF shown by Eddington (1916) and many different anisotropic DFs obtained by many other outstanding 
astronomers (e.g., Camm 1952; Bouvier 1962, 1963; Kuzmin & Veltmann 1967, 1973; Veltmann 1961, 1965, 1979, 1981; Kent 
&Gunn 1982; Dejonghe 1986, 1987; Dejonghe & Merritt 1988). Recently, a method was presented by Jiang & Ossipkov (2007) 
for finding anisotropic distribution functions for spherical galaxies. This is an combination of Eddington's (1916) formula and 
Fricke's (1952) expansion idea. Of course, they can be also regarded as simply an extension of the idea of Eddington. 

Fricke's (1952) expansion idea is that DFs which are products of the two powers of the energy and the square of the 
angular momentum about the axis of symmetry correspond to densities which are proportional to products of the potential 
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and the radial coordinate for axisymmetric systems. Hence the DF for the system can be obtained by first expressing the 
density as a function of the potential and the radial coordinate, and then expanding as a power series. According to the 
maximum entropy principle, there are an infinity of the most probable two-integral DFs for a given mass distribution (see 
Dejonghe 1986). Thus there may be an infinity of two-integral DFs corresponding to any given mass density in axisymmetric 
stellar systems. 

In the literature, there are a number of different axisymmetric models for galaxies (e.g., Miyamoto 1971; Bagin 1972; 
Miyamoto & Nagai 1975; Nagai & Miyamoto 1976; Kutuzov & Ossipkov 1980, 1986, 1988; Evans 1993, 1994; Kutuzov 1995; 
Jiang 2000; Jiang & Moss 2002; Jiang et. al 2002; Jiang & Ossipkov 2006; Ossipkov & Jiang 2007) and different integral 
transformation techniques are used to obtain the solution of the problem of findng the two-integral DFs of axisymmetric 
stellar systems (e.g. Lynden-Bell 1962; Hunter 1975; Kalnajs 1976; Dejonghe 1986) but there is a fatal obstacle of requiring 
not only the validity of these transformations of the density but also the complex analyticity of a density-related integral 
kernel to complex arguments. It is a great progress that the contour integral of Hunter & Qian (1993) is used to find the 
two-integral DFs for axisymmetric systems. However, this is still a complex integral suitable for analytic densities and its 
singularities are required to satisfy some conditions. 

Therefore the above ideas of Eddington and Fricke are still very useful to the problem of finding the DFs. The main 
aim of this paper is to present some formulae for finding two-integral DFs in axisymmetric systems by an combination of 
both the Eddington formula and the Fricke expansion (see above). The fundamental integral equations of the problem of 
finding the two-integral DFs are given in Section 2 and some new formulae of the two-integral DFs for stellar systems with 
known axisymmetric densities are shown in Section 3. These DFs depend only on the two integrals of the energy and the 
magnitude of the angular momentum about the symmetry axis. These formulae are also an extension of those shown by Jiang 
& Ossipkov (2007) for finding anisotropic distribution functions for spherical galaxies. A type of two-integral DF which is a 
sum of products of functions only of the energy and powers of the magnitude of the angular momentum with respect to the 
axis of symmetry is derived in Section 3.1 and another, which is a sum of products of functions only of a special variable 
and powers of the magnitude of the angular momentum about the axis of symmetry, in Section 3.2. More general formulae 
are given in the last part of Section 3. Various formulae of the velocity dispersions for such models of these DFs are also 
shown in all the three parts of Section 3. For the gravitational potentials having no upper bound and tending to +oo, the 
similar formulae of the DFs are also shown in Section 4. Their application is given in Section 5. Section 6 is a summary and 
conclusion. 



2 THE FUNDAMENTAL INTEGRAL EQUATIONS 

Assume that $ and E are, respectively, the gravitational potential and the energy of a star in a stellar system. As in BT, 
one can choose a constant $o such that the system has only stars of the energy E < $o, and then define ip = — $ + $o 
and £ = — E + 3>o. In this system, the two physical quantities tp and s are usually called the relative potential and energy, 
respectively. Obviously, e = is a relative energy of escape from the system. Given a stellar system with the relative potential 
tp — ip(r), its mass density p — p(r) can be obtained by using Poisson's equation and its distribution function / = /(r, v) 
satisfies the following integral equation 



where r is a position vector, v is a velocity vector, G is the gravitational constant. The cylindrical polar coordinates (R, tp, z) 
are generally used with the z-axis being that of symmetry for an axisymmetric system. The velocity in cylindrical polar 
coordinates (R,ip,z) is usually denoted by v = (vr,v^,v z ) and L z is the component of angular momentum about the 2-axis. 
Then L z = Rv$ and it is well known that the relative energy e and the 2-axis angular momentum L z are two isolating 
integrals for any orbit in the axisymmetric system. Hence, by the Jeans theorem, the DF of a steady-state stellar system in 
an axisymmetric potential can be expressed as a non-negative function of e and L z , denoted by f(e,L z ), and then for an 
axisymmetric system, (1) can be rewritten as 



Let /+(e, L z ) = [/(e, L z ) + f(e, —L z )\/2. Then, by e = tp — (1/2)(vr + + v z ), the integral given by (2) can be expressed as 



since the system has only stars with e > 0, that is, f(e,L z ) = for £ > 0. This implies that a given density determines 
/+(e, L z ) which is just the part of the DF that is even in L z . Hence /+(e, L z ) is usually called the even DF. 




(1) 




(2) 




(3) 
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Once f+(s, L z ) is known, f(e, L z ) can be obtained under some suitable assumptions such as the isotropy of the two- integral 
DFs (see BT) and the maximum entropy of the most probable two-integral DFs (see Dejonghe 1986), and further p can be 
easily calculated by integration and tp by solving Poisson's equation for the axisymmetric system. The inverse problem that 
is now investigated is how to derive the two-integal DF f+(e,L z ) from the density p for any axisymmetric system. Different 
classes of the two-integal DFs will be below shown, which are derived from axisymmetric density profiles for galaxies, by 
combining some functions only of e (or Q = e — L 2 /(2R 2 )) with some functions of the form L 2nl3 " where R a is a scaling 
radius, n is an integer greater than —2 and /?„ is a constant such that n/3 n > — 1. 



3 TWO-INTEGRAL DFS FOR AXISYMMETRIC SYSTEMS 

In this section various formulae for the even DFs of axisymmetric systems are obtained from axisymmetric density profiles of 
different forms and their radial and rotational velocity dispersions are expressed in a simple form. 



3.1 DFs of the form £ L 2 z n h n (e) 

71 = 

Note that the integral on the right side of (3) is in fact a function of the relative potential ip and the radial coordinate R. 
Hence, by (3), the mass density p can be regarded as a function depending on the relative potential ip and the radial coordinate 

m m 

R. Let p = p(r) be below denoted by p{ip,R). Assume that f + (e,L z ) — ^2 L 2n h n {e) and that p(ip,R) = R 2n p„(ip). Then, 

71=0 71=0 

by (3), it is easy to see that 

m m /^o« + l/2p2n /■</> 

E fl2 "j?"W = E 2 n + l / h n (e)(ip~er + ^de. (4) 

n=0 n=0 ^° 

It follows from (4) that 
->n+l/2 ri> 



PnW = 2n + 1 I h n {e){i> - eT +1/2 ds. 



(5) 



Assume that (d? p~ n (i>) /dtp J )^ = o = for all j € {0,1, ... , n} and all n G {0, 1, 2, ... , m}. By taking the (n + l)st derivative of 
(5) and using Abel's integral equation, one can get 

n( ) (27r)3/22"F(n + i) de J d^ ^/T^P ( > 

for n = 0, 1, 2, ... , m. Hence it can be easily known from (6) that 

i ^ L 2 z n d r d- +1 p n w d^ 

J+( ' z) (27r)3/ 2 ^2"r(n+i)de7 d^+i [ > 

m 

for e > 0, corresponding to the axisymmetric density of the form p(tp,R) = ^ Pn{"4>)R 2n ■ Also, (7) can be equivalently 

n=0 

rewritten as 



1 f d n + 2 p n W dj> 1 ^rf" + 1 pn(V) 

J+l£ ' zJ (2tt)3/ 2 4^ 2™r(n+i) 7 <#™+ 2 v^ 3 ? 



Ve - i> y/e \ dip 



•0=0 



(8) 



Furthermore, if it is assumed that (d 3 p n (tp) /dip 3 )^, =0 — for all j G {0, 1, . . . , n + 1} and all n G {0, 1, 2, . . . , ra}, then, for 
e > 0, (7) can be expressed as 



1 ."\ jj 2n f e 



It is worth mentioning that Evans (1994) used Abel transforms to get a similar formula of the above even DF in the case 
for m = 1. His formula of the even DF has just two previous terms and is formally expressed as 

r2 d T dh{4>) dip 



j_d_ r rfpow dxp l 2 z d r 



(10) 



dxp (e-i/>) 1/2 2V8n2deJ Q d^ (e - V) 3 / 2 ' 

The real integral in the second term on the right side of (10) is divergent for the stellar systems such as the Lynden-Bell 
model considered below. 

By (7), the velocity dispersions a 2 R (ip,R) and a 2 (ip,R) can be also found as follows: 
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and 

m 

4(^i?) = ^ r - S y^(2n + l) J R 2 " / (12) 



n=0 



for any DF derived from the axisymmetric density of the form p(ip,R) = ^2 Pn(4>)R 2n , here and below, v v represents the 

n=0 

mean rotational velocity which can be calculated under some suitable assumptions. It can be also known that these dispersions 
(11) and (12) can be obtained directly according to Hunter's (1977) formulae as follows: 

a%(iP,R) = a 2 (^,R) = tf^]{jf o PW',W (13) 
and 

R) = pJ^R) J, OR *+ ~ v "- (14) 



3.2 DFs of the form £ L 2 "g n (Q) 

n=0 

A more general expression for the integral in the right side of (1) can also be derived. Put Q = e — L z /(2R a ), where R a is a 
scaling radius, and assume that the DF is dependent on Q and L z , denoted by f{Q, L z ), and that the system has only stars 
with Q > 0, or equivalently, / = for Q < 0. Obviously, Q — > e as R a — > oo. Then, for an axisymmetric system, (1) can be 
expressed as 

s>2„ 



-ira(*£)-&-«'-W /< * I -*" v - 

By changing the variables of the integral in (15), it follows that 

r * f /•v /2 (*-<3)/( 1 +- R2 /fla) 



/"A 

■ = 4tt / / 



f+(Q,L z )dv4, 



dQ, 



(15) 



(16) 



where f+(Q,L z ) = [f(Q,L z ) + f(Q,-L z )}/2. Naturally, f+{Q,L z ) is the even part of f{Q,L z ). Suppose that f+(Q,L z ) = 

m 

^2 L 2n g n (Q)i and that the mass density has the following form 



pty, R) = Y, MV0^ 2 7(i + R 2 /Rlr +1/2 - 

Then (16) can be expressed as 

R 2n Mip) ^ ^2 n+1 ' 2 R 2n '+ 



E47TZ 
. . .. (2n + l)(l + R 2 /R 2 )"+ 1 / 2 

n — n — 

It is easy to see that equation (18) gives 

47r2 n + l/2 ri> 



Jo 



94Q)^-Q) n+1/2 dQ 



2n + l 



g4Q)^~QT +1/2 c 



(17) 



(18) 



(19) 



Assume that (d 3 p n (tp) /dtp 3 )^=o = for all j € {0,1, ... , n} and all n £ {0, 1, 2, ... , m}. By taking the (n + l)st derivative of 
(19) and using Abel's integral equation, one can obtain 



9n(Q) 



1 



d 



1 p„{i>) dip 



(20) 



(27r)»/ 2 2"r(n + i) dQj Q VQ^P 

for n = 0, 1, 2, ... , m. Hence, corresponding to the axisymmetric density p(ip, R) defined by (17), it can be readily shown from 
(20) that the DF is obtained as 

{2Tv)- 3 ' 2 L 2 z n d f Q d" +1 p n (i,) dip 



f+(Q,L z ) = J2 



2«r(n+i) 



dip n+1 VQ~^ 



(21) 
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or equivalently, 

" Q d"+ 2 p n (V) dxp 1 (d n+1 p n {i>) 



m /o-1-3/2r2n r 



i/>=0 



(22) 



for Q > 0. Furthermore, if (cP p„(^)/dip j )^ = for j € {0, 1, . . . , n + 1} and all n £ {0, 1,2, ... , m}, then for Q > 0, (21) 
can be rewritten as 

{ ] ~h ^ n nn+\) I VQ -^- (23) 

Of course, (21) and (23) coincide with (7) and (9), respectively. In other words, (7) and (9) are, respectively, limits of 
(21) and (23) when R a — > co. 

Dejonghe (1986) studied the mass densities separable in ip and R and gave the following formula of the even two-integral 
DFs: 

/3 +ioo 

r( P +l) £ ^/ 2 1 f M R ^{g}[L 2 /(2e)]-^ 2 
U(£ ' z) ~ 23/^77 27 ri J r(l/2-/3/2)/3(p-l/2+/3/2)^ l24j 

/3 -ioc 

for any given mass densities p(tp,R) = tp p g(R), where i = v 7 — 1, p > 3/2, f3o is a suitable constant, ^ffl R ^/3{g} represents the 
Mellin transformation of the function g(R). Let a and b be two constants. Assume that —a — b is not a natural number. When 
g(R) = i? 2 7(l + i? 2 ) a+f \ (24) is written by Dejonghe (1986) as 

F(v+ lk p ~ 3/2 1 1 L 2 

M°> L ')= Xrt + b) ^ b >>-2>2>& (25) 

where H(a, 6, c, d; x) is defined by 

with the contour C such that —a are on its left side and b on its right side. In the case that a + d and b + c are not negative 
integers, the complex integral (26) can be calculated and expressed as follows. When < x < 1, if a — cisa nonnegative integer, 
then H(a, b, c, d; x) = 0, or else H(a, b, c,d;x) = x a 2 Fi(a + b, 1 + a — c; a + d; x)F(a + b)/[r(c—a)T(a+d)]; when x > 1, if 6 — d is a 
nonnegative integer, then H(a, b, c, d; x) = 0, or else H(a, 6, c, d; x) = x~ 6 2 fi (a+6, 1 + 6 — d; 6 + c; l/x)r(a + fc)/[r(d — fc)r(fe + c)]. 
Here, 2 Fi is a hyergeometric function. In particular, H(n, 1/2, p — 1/2, 1/2; x) = x" 2 F 1 (n + 1/2, n — p + 3/2; n + 1/2; x)/F{p — 
n - 1/2) for < x < 1, and H(n, 1/2, p - 1/2, 1/2; x) = for x > 1. Therefore, when a = n and & = 1/2, the DF given by (25) 
can be rechanged as 

t , r ^ r(p+l)^- 3 / 2 (LlY 1 3 1 L 2 

^ £ ') = 2»/^(p-n-l/2) 2Fl(n+ 2^-^ + 2 ;n+ 2'2F ) 

r(p + l)e p - 3/2 (hi V g r(n - p + fc + 3/2) f L 2 x 



23/2 7r r(p-7i- 1/2) ^26/ r(n-p + 3/2)fc! ^ 2e 

7r- 1 r(p+ l)L 2n e p - n - 3 / 2 f 1 _£z_ 



2 n+3/2 r ( n + l)r(p _ n _ 1) ^ 2e 



^T( P +1)L^ ( £2\ P-n-3/2 



2«+3/2 r ( n+ i)r(p- n - i) 



(27) 



for e > L 2 /2 and f+{e,L z ) = for e < i 2 /2. It can be easily found that the DF (27) is the same as obtained by use of (23) 
when p — n > 1. 

Similar to those in Sect. 3.1, the velocity dispersions a R (ip,R) and a 2 (ip,R) can be also found to be of the following 
forms 



*ti> — ^77 Y'. — ,>.",^+ 1/2 J Mv')dv' (28) 



R 2 

P {i>,R) ^ (1 + R 2 /R 2 Y 



n=0 



and 



1 m f? 2 ™ 

R ^pJ^R) + (1 + flV ^)^ I WW-^ (29) 

n— 

for any DF derived from the axisymmetric density p{tp,R) defined by (17). 
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3.3 Miscellaneous DFs 



One can also obtain more general formulae than (7) and (21). Assume that Q = max(Q,0). Then it can be further shown 
that the DFs of the form 



f+(e,Q,L z )= J^BmL 



2nf) lr , 



re d a i" + 1 p„(» dij, , i_ ( d^±nX±)\ 

JO di/> a ln + 1 (e-V) 01 " \ di, a ln J 



+ E B 2n L 2n ^« 

n=0 



JO dV<"2„ + l (Q_V) Q 2n Q°2„ ^ dV-"2„ 



(30) 



correspond to an axisymmetric density of the form 

m m 
n=0 n=0 

with n/3i„ > —1, where Bin = [(27r) 3/,2 2 n/3in r(7i/?m + 1/2)T(1 — ai„)] _1 , a™ = n(3i„ — cn„ + 3/2 and a,i„ is a non-negative 
integer such that < cn n < 1 for i = 1, 2 and n = 0, 1, • • • , m. 



Finally, the velocity dispersions a R (ip,R) and a v (ip,R) can be also obtained as 



^ = ^E ^/ p4W + 7TT 



I p n {ip')dip' 
Jo 



and 



R 2 /Rl) n ^n+l/2 

(2n[3 2n + l)R 2n ^ 



R 2 /R 2 a) n ^+3/2 



I p n (tp')dip' 
Jo 



(32) 



(33) 



for any DF derived from the axisymmetric density p(tjj,R) given by (31). 



4 MODELS WITH GRAVITATIONAL POTENTIALS HAVING NO UPPER BOUND 



The axisymmetric gravitational potential <E> now has no upper bound and tends to +oo at large distances from which escape 
is impossible. Thus one usually denote by f(E, L z ) the two-integral DF of a steady-state stellar system with the axisymmetric 
potential. Suppose that the system has only stars of E > 0. Then, by using the even two-integral DF f+(E, L z ) = [f(E, L z ) + 
f(E, —L z )]/2, the fundamental integral equation (1) can be changed as follows: 



■=~r\[ 

R i<f Jo 



R^/2(E-<S>) 



f+(E,L z )dL z 



dE. 



As in Section 3.1, one can obtain the similar even DF 



{ _ irL 2n 



T(n+i) 



— r 

dEj E 



d n+1 p~„{®) d$ 



(34) 



(35) 



for E > 0, corresponding to the axisymmetric density of the form p($, R) = Y] p n ($)J? 2n under the assumption that 

n=0 



S^ + oo dQ3 E^ + oc 



d j p n ($>) d$ 



d$J y$ _ E 







(36) 



for all j G {0,1,..., n} and all n £ {0, 1, 2, . . . , m}. Furthermore, if it is assumed that the condition (36) holds for all 
j € {0, 1, . . . , n + 1} and all n £ {0, 1,2,..., m}, then, for E > 0, (35) can be expressed as 



f+(E, L z 



^ m 
(2^72 51 



(-1)"L 2 
2™r(n + 



n 7- 2n 



d n+2 p n {$) rf$ 

d$™+ 2 E' 



(37) 



It can be further shown that (35) and (37) are at least formally in accordance with the contour integrals given by Hunter and 
Qian (1993) for the gravitational potential tending to +oo at large distances. 

Put Q = E + L 2 Z /(2R 2 ) for the system with gravitational potentials having no upper bound. Similar to those given in 
Section 3.2, one can also obtain the even DF 



f+(Q,L z )=J2 



(2tv)- 3 / 2 L 2 z " d 
(_2)»r(n+i)dQ 



d n+1 p n {$) d$ 



(38) 
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for Q > 0, corresponding to the axisymmetric density R) of the form defined by 

m 

R) = J2 /Sn(<£>)i? 2 7(l + R 2 /Rl) n+1/2 (39) 



n=0 

under the assumption that 

hm — y = and lim Q / — ' , = (40) 

hold for all j £ {0, 1, . . . ,n} and all n € {0, 1, 2, . . . ,m}. Furthermore, if one assume that the condition (40) holds for all 
j € {0, 1, . . . , n + 1} and all n £ {0, 1, 2, ... , m}, then for Q > 0, (38) can be rewritten as 

/+(0 ' Lz) "i(-2)"r(n + i)y Q V *-Q- (41) 

As in Section 3.3, one can further show the more general DFs of the form 

m 

f + (E,Q,L z )= E (-ir^B^L^^J^ d £ a M*> 

n=0 
m 

+ E (-i)« 2 »-i BanL ;»fc.^jj»^^_ (42) 

n= 

corresponding to an axisymmetric density of the form given by 

m m 

p($,J?) = ^p„($)i? 2 "^+^ / 3„($)i? 2 " fe V(l + ^ 2 /-Ra) n ' 32 ' l+1/2 I (43) 

n=0 n=0 

where Q = max(Q,0), a in , B in , Q 4n and (3 irl are the same as in (30) for i = 1,2. 

By (34), it can be also found that, in the system with only stars of E > 0, the even DF 

/+(£, L.) = L 2 " +1 cxp ^£ - g|) (44) 

corresponds to an axisymmetric density of the form 

47^(2n)!! J Rg ( " +1) i^ 2 " +1 e-"^ I ' 
a(R%a + /3Ji 2 )"+ 1 

for any gravitational potential tending to +oo, where a and /3 are nonnegative constants, 7?o is a positive constant, n is a 
nonnegative integer, and (2n)H = 1 • 2 • 4 • • • (2n — 2) • (2n) when n is a natural number and (0)!! is defined to be equal to one. 
It is very remarkable that (44) can be recovered from the complex contour integral given by Hunter and Qian (1993). It can 
be below used to find the odd part of the DF that corresponds to some assumed rotational velocity {v<t>) for the axisymmetric 
stellar systems. 



rf*. R) = " XZ , gfk^ («) 



5 APPLICATION TO THE AXISYMMETRIC CASES 

Binney's (BT) logarithmic model has infinite mass and its gravitational potential is of the form 

$(R,z) = ±v 2 ln(l + R 2 + ^j , (46) 

where vo is the constant circular velocity in the equatorial plane at large distances, q is the axial ratio of the spheroidal 
equipotentials. Obviously, this gravitational potential has no upper bound. The density derived from (46) is 

p(R,z) = {2[(1 - g 2 )i? 2 + l]e- 4 */"o + {2q * _ 1)c -»/«g} . (47) 

Then, by (37), one can obtain the even DF corresponding to (47) as follows: 

ME, L t ) = {2 9 / 2 [(l - q >)L\ + <l^vl\e-™l< + (2g 2 - l^e" 2 ^ } , (48) 

which is the same as found by Evans (1993) using Lynden-Bell's (1962) method. This kind of solution was known earlier to 
Toomre (1982) and published first by Miller (1982). It can be easily found that the mass density (47) of the Binney model 
is positive in the position space only when the axial ratio q is not less than l/y/2. It is easy to see that, in the prolate case 
when the axial ratio q is greater than 1, the DFs (48) must be negative at some points in the physical domain and so they 



© 2007 RAS, MNRAS 379, 1133-?? 



1140 Z. Jiang and L. Ossipkov 



4.0 — r 

3.5 — 
3.0 — 
2.5 — 
kj 2.0 — 
1 .5 — 
1 .O — 
0.5 — 



(b) 



O.O -f 



-^1 



4.0- 
3.5- 
3.0- 
2.5- 
2.0- 
1 .5- 
1 .O- 
0.5- 

o.o— f- 



(a) 



O.O 0.5 



2.0 2 

E 

4.0 — i- 
3.5 — 
3.0 — 
2.5 — 
2.0 — 
1 .5 — 
1 .O — 
0.5 — 



(c) 



C3 



O.O — f- 



-^1 



O.O 0.5 



2.0 

E 



3.5 4.0 



O.O 0.5 



2.0 

E 



Figure 1. The contours of the DFs given by (48) with three different axial ratios when vq is set to be one. Subfigures (a), (b) and (c) are 
for q = 1, q = 0.9 and q = 0.8, respectively. The thin solid curves are isocontours and the dotted curve is the boundary of the physical 
domain. Successive levels differ by factors of 0.4. 



are not the real two-integral DFs for the stellar system. When the axial ratio q is equal to one, the model is spherical. The 
model is flattened only when 1 > q > l/\/2 w 0.707107. Figure 1 illustrates the contours of the DFs given by (48) with three 
different axial ratios. One of them is for the spherical case and the other two for the flattened one. 

Note that if p and f + (E,L z ) in (34) are replaced by pR{v^,} and L z f_(E,L z ), respectively, then (34) becomes an equation 
for L z f_(E,L z ), that is, 



47T 



L z f_{E,L z )dL z 



dE, 



(49) 



where f_(E,L z ) is usually called the odd DF given by f_(E,L z ) — [f(E,L z ) — f(E,—L z )]/2 for the stellar systems. This 
property was first found by Lynden-Bell (1962) and then applied by Hunter and Qian (1993) into calculating the odd part of 
the DF for the Binney model under the assumption of the rotational velocity {v$) having the rotation law: 



M =v t R 2 /(R 2 t +R 2 ), 



(50) 



where and R* are constant velocity and length scales, respectively. (50) was one of the rotation laws considered by Evans 
(1993). Using the Hunter and Qian's (1993) contour integral formulae of the even DF, Hunter & Qian (1993) first derived a 
contour integral of the odd DF from the rotational velocity {v$) for the Binney model and this idea was then extended by 
Jiang (2000) to the odd DF for a more general model. By combining (44) and (49), it can be also found that, in the system 
with only stars of E > 0, the odd DF 



f_(E, L z ) = sgn(L z )L 2 z n exp -aE 



2R 2 

corresponds to an axisymmetric density p($, R) of the form 
47r(2n)!!flg (n+1) fl 2 "e- Q ' 1, 



a(R 2 a + /3i? 2 )"+! 



(51) 



(52) 



for any system with a gravitational potential tending to +oo, where a, f3, Ro, n and (2n)!! are the same as in (44). If one 
assumes that the rotational velocity (t^) of the Binney model satisfies the rotation law 

= v*R 2(n+1) /{Rl + R 2 ) n+1 (53) 
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for any nonnegative integer n, where v* and J?* are the same as in (50), then, combining (47), (51) and (52), one can find 
that the odd DF is given by 

u„sgn(Lz) / 2 L 2 Z ( 4E\ 2 2 / 4E 



f-(E, L z ) = ^J G y v 2 \ 16(1 - Ott exp ( --y ] + 8[1 - (n + !)#(! - 9 )] exp 



o 



2 2j / L z \ 2j / 4£ 2L 2 \ ,„ 2 / 2£ 

-sEMnl^) exp^-^j+^-Dexp^--^ 

7=0 



n/2 2fc-l „ ■ „ ■ / 

2 2j / L z \ 2j / 4£ 2L 2 



+ 16^(l^ 2 )sgnH5:5: 7 i- -g- exp 



(2j)i! "V «s 

fc=l 3=0 
3=0 

j=0 ) 



for any nonnegative even number n and that the odd DF is expressed as 
««sgn(L z ) 

2 2j / L z \ 2j ( 4E 2L 2 \ , ,„ 2 ^ / 2_B 



^) = ^TOT 1 16(1 - M exp (-^) + 8[1 - (n + 1)^(1 - ^)] exp (-*§ 



E2 iJ ( L z Y J 4E 2Lt \ ,„ 2 ,x / 2S 

T^TTf T^" ex P 2" - T^TT + 2 ? ~ 1 ex P 2" 

(2j)!! \R*v J \ vl Riv^l \ v o 

j=0 

(n-l)/2 2fc „. „. / 

+i6^(i- 9) x: E(2^(^j -p^-^-^i 

fe = j = o w ' \ u u . 

2 N ( "v^ /2 2 4 ^+ 2 / L z \«+ 2 / 4E 2Ll " 
i=o 

3=0 ^ ^ 

for any positive odd number n. This also means that there are an infinity of the DFs for any given axisymmetric stellar 
potential. (53) is obviously an extension of (50), and when n = 0, (54) is the same as given by Hunter and Qian (1993). 
The well-known Lynden-Bell (1962) model has finite mass and its relative potential and density are given by 

4>(R, 2) = [(R 2 + z 2 + l) 2 + aR 2 ]- 1/4 , (56) 



(57) 



(3 + a) - 5a (l + |) R 2 i> 
where a is a flattening parameter. By using (7) and (57), it follows that the even DF is given by 



2 3 / 2 tt 2 \ 143 

which is in fact as the method of Fricke (1952). It is here necessary to explain the different definitions of the gravitational 
potential. In the paper written by Lynden-Bell (1962), (56) is called the gravitational potential of the Lynden-Bell model. 
Due to the use of the concept of the relative potential, the gravitational potential defined by Binney and Tremaine (BT) in 
fact differs by a factor —1 from that given by Lynden-Bell (1962). The Binney and Tremaine definition of the gravitational 
potential is used throughout this paper and so (56) is a relative potential. 



6 CONCLUSIONS 

Few galaxies are even nearly spherical. Thus it is a natural idea to explore some important properties of real galaxies by 
employing the cylindrical polar coordinate system (R, ip, z) with the center on the galactic nucleus and the z-axis being that 
of symmetry of the galaxy. However, at least two involved factors require being mentioned as follows. One is that some possible 
orbits in many real galaxies can be easily described by studying a two-dimensional problem. With the help of the conservation 
of the angular momentum about the symmetry z-axis, this problem can be directly reduced from the analysis of the orbits 
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in the three-dimensional space occupied by axisymmetric galaxies. The other is that, on the analogy of anisotropic DFs for 
spherical systems, the two-integral DFs for some axisymmetric systems can be also found to model the typical behaviours of 
the dynamical quantities of galaxies considered. 

Some formulae of the two-integral DFs can be obtained for stellar systems with known axisymmetric density as a sum of 
products of functions only of the potential and a special function (or power) only of the radial coordinate, i.e. these DFs are 
a sum of products of functions only of a special variable (or the energy) and a power only of the magnitude of the angular 
momentum about the axis of symmetry. They come from an combination of the ideas of Eddington and Fricke and they are 
also an extension of those shown by Jiang and Ossipkov (2007) for finding anisotropic distribution functions for spherical 
galaxies. As an analogue for spherical models, the product of the density and its radial velocity dispersion can be also expressed 
as a sum of products of functions of the potential and of the radial coordinate. But the expression of its rotational velocity 
dispersion formally differs from that of its radial velocity dispersion. It can be further found that the density multiplied by 
the difference between the dispersion of its rotational velocity and the square of its mean rotational velocity is equal to a sum 
of products of functions of the potential and of the radial coordinate. The similar formulae of the two-integral DFs for the 
gravitational potentials without upper bound are as well in accordance with the complex contour integral ones given by Hunter 
& Qian (1993). These expressions for axisymmetric systems can be used to obtain the even DF of Binney's (BT) logarithmic 
potential although Evans (1993) derived it using Lynden-BelPs (1962) method. An infinity of the odd DFs for the Binney 
model can be also found under the assumption of the laws of the rotational velocity. For the well-known Lynden-Bell (1962) 
model, these analogues degenerate into the method of Fricke (1952). It is worth mentioning that such analytic procedure 
to determine the DFs can be also applied to the prolate Jaffe models given by Jiang & Moss (2002) for a good numerical 
approximation of the two-integral DFs for the stellar systems. 

One can finally know that it is a shortcoming of all the two-integral models that the radial velocity dispersion is equal 
to the vertical velocity dispersion. This is because it is well-known that in real axisymmetric stellar systems, the velocity 
dispersion in the radial direction is not equal to the velocity dispersion in the vertical direction, meaning that the DFs of 
the real systems must actually depend on three integrals of the motion (one of them being non-analytic in general) rather 
than two (of course, there is then no unique solution for the even part of the DF). To overcome this shortcoming of the 
two-integral models, some extensions of two-integral DFs have been studied to construct three-integral DFs for particular 
orbital families in flattened axisymmetric systems (Evans, Hafner & de Zeeuw 1997) and for separable axisymmetric Stackel 
potentials (Famaey, Van Caelenberg & Dejonghe 2002). 
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